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Plectoneme creation reduces the rotational friction of a polymer 
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Abstract. - The torsional dynamics of a semiflexible polymer with a contour length L larger than 
its persistence length L p that is rotated at fixed frequency ujq at one end is studied by scaling 
arguments and hydrodynamic simulations. We find a non-equilibrium transition at a critical 
frequency to*: In the linear regime, uio < <^«, axial spinning is the dominant dissipation mode. 
In the non- linear regime, uio > <^», the twist-dissipation mode involves the continuous creation of 
plectonemes close to the driven end and the rotational friction is substantially reduced. 
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Introduction. — The properties of semiflexible poly- 
mers have received continuous interest and attention from 
biophysicists for the last two decades. This is so since 
they constitute a minimal yet realistic mechanical model 
for biological macromolecules such as DNA and filamen- 
tous proteins. The framework of statistical mechanics pro- 
vides various equilibrium properties of semiflexible poly- 
mers such as the force-extension relationship [1], which 
have been successfully tested through detailed comparison 
with experimental data obtained from newly developed 
micromanipulation techniques employing optical tweezers 
and magnetic beads [2] . The importance of torsional elas- 
ticity has been recognized in the context of DNA supercoil- 
ing [3] . Driven by experiments where both stretching and 
twist of single biopolymers could be controlled [4], theo- 
ries were generalized to include twisting effects as well and 
provide quantitative agreement on the static level [5]. 



The in-vivo functioning of DNA involves non- 
equilibrium dynamic twist effects that mostly have to do 
with the activity of various DNA- processing proteins [6] . 
For example, in replication, the process by which a DNA 
chain is duplicated into two identical daughter strands, 
the helical nature of DNA requires daughter strands to 
unwind and thus the mother strand to rotate around its 
axis. This large scale motion has been considered a con- 
ceptual obstacle; it was resolved by a simple calculation 
demonstrating that the friction dissipation associated with 
simple axial-spinning of DNA (like in a speedometer ca- 
ble) is rather small compared to typical biological free en- 
ergies [7]. In transcription, the process by which the DNA 



informational content is copied into a continuously grow- 
ing RNA chain, it was suggested that a long nascent RNA 
chain might (either due to its own hydrodynamic friction 
or via anchoring to some other cellular component) pro- 
vide enough rotational resistance to force the DNA strand 
to rotate around its own axis [8] . Experimental in-vivo and 
in-vitro studies found large degrees of DNA supercoiling 
upon transcription, being positive in front and negative 
behind the transcriptional complex [8] . This finding, how- 
ever, is at odds with the above-mentioned simple axial- 
spinning scenario [7] , since the rotational friction of bare 
DNA seems not large enough to induce sufficient torsional 
stress needed to induce supercoiling [9]. In an effort to 
reconcile the conflicting experimental results, Nelson in- 
troduced the notion of static bends along DNA, which 
substantially increase the rotational friction [9]. So far, 
the issue of the rotational friction of a semiflexible chain 
is still a matter of interest and debate. Single-molecule 
experiments on rotating DNA failed to observe the static- 
bend-induccd anomaly of the rotational friction, but this 
could possibly be related to the high rotational frequency 
employed in the experiments [10]. Several studies did not 
provide a clear-cut answer as to what the influence of static 
DNA bends is and whether other factors present in in-vivo 
studies are needed to induce substantial supercoiling of 
DNA under twist injection [11-14]. 

In this Letter we use non-equilibrium scaling arguments 
and hydrodynamic simulations and consider a fundamen- 
tal scenario, namely a homogeneous semiflexible polymer 
in the stationary limit, that is axially rotated at one end 
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at frequency ujq with the other end free, as schematically 
depicted in Fig. 1. In the absence of shape fluctuations, 
i.e. at zero temperature, an elasto-hydrodynamic insta- 
bility occurs at a critical driving frequency lj c that sepa- 
rates a twirling from a whirling regime [15]; for luq < lo c 
the rod stays straight and undergoes simple axial spinning 
(twirling), while for loq > u> c the rod buckles and displays 
a combination of axial spinning and rigid-body rotation 
(whirling). Including thermal effects, elastic shape fluc- 
tuations of the rod round off the transition and shift the 
instability to lower frequency [16] . In the present study, we 
consider the opposite limit where the chain contour length 
is much larger than bend and twist persistence lengths and 
the twirling-whirling instability is completely washed out 
by thermal fluctuations. This is relevant to transcription- 
driven DNA supcrcoiling, but has not been studied as an 
elasto-hydrodynamic problem for a homogenous scmiflcx- 
iblc chain before. As our main result, we find in addition 
to the well-known axial-spinning regime, realized for low 
driving frequency ujq < u>* , a novel plectoneme dominated 
regime at high frequencies luq > cj*. In this regime, twist is 
converted locally into writhe (in the form of plectonemes) 
close to the driven end and then diffuses out to the free 
end. For sufficiently long chains, the crossover frequency 
c<j» is much larger than the twirling-whirling threshold io c . 
Quite surprisingly, in the plectoneme regime the overall ro- 
tational friction is significantly reduced as compared to the 
axial-spinning scenario, i.e., the plectoneme dissipation 
channel transforms the injected twist into writhing motion 
at very low frictional cost. This mechanism could be bi- 
ologically relevant since it shows how transcription might 
produce positively super-coiled DNA structures that could 
favorably interact with negatively super-coiled nucleoso- 
mal structure ahead of the RNA polymerase [6] even in 
topologically open systems and at very small energy ex- 
penditure. 

Model. — Within linear elasticity theory, the elastic 
energy of an isotropic rod of contour length L is given by 



ds 



(1) 



where A and C are the bend and twist rigidities, k(s) = 
\d 2 r\ is the local curvature with r(s) the position of the 
rod centerline and d s denoting the derivative with respect 
to the arclength s. The local twist density fi(s) is defined 
by Q(s) = d s <t>{s) where <f> is the rotational angle about 
the local tangent d s r. Note that in cq. (1) includes 
also the twist contribution due to geometric torsion, i.e. 
changes of fa[s) are also generated e.g. by out-of-planc 
bending [17,18]. 

In our simulations, a filament is modelled as a chain of 
N + 1 connected spheres of diameter a. The total en- 
ergy of the system, E = E st + E e \ + Elj, includes a 
stretching contribution ensuring connectivity of spheres, 
E st = K/2J2j=i(\ r j+i ~ r j\ _ a ) 2 i where Tj is the po- 
sition vector of monomer j, and a truncated Lennard- 
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Fig. 1: A semiflexible polymer is rotated at frequency ljq at a 
fixed end and exhibits rotation at frequency ujl at its free end. 
The three twist-dissipation channels are (A) axial-spinning, (B) 
solid-body rotation, and (C) plectoneme-creation. 



Jones potential to account for chain self- avoidance, where 



E L j = £LjE 4 <j[( a / 



< a. 



12 — 2(a/ry) 6 ] applied only for 
The local elastic translational 
force Fj and torque about the local tangent Tj act- 
ing on each sphere are calculated using the variational 
method described previously [18], leading to the coupled 
Langevin equations: dtYi = Y^j^i A**?' ' + an d 
dt4>i = (irrja 3 ) Tj + Sj(t), where fa is the rotation an- 
gle around the bond vector r^+i — r.;. Hydrodynamic 
interactions between two spheres i and j are included 
via the Rotne-Prager mobility tensor /j,^ [19]. For the 
translational self-mobility of the spherical monomers with 
diameter a we use the standard Stokes formula flu = 
1 / (37r?7a) = /iol- The vectorial random displacements £(£) 
and E(£) model the coupling to a heat bath and obey the 
fluctuation-dissipation relations implemented numerically 
by a Cholesky factorization [19]. For the numerical inte- 
grations we discretize the Langevin equations with a time 
step A and rescale all lengths, times and energies, leading 
to the dimensionlcss parameter A = Afc^T/in/a 2 . We set 
the twist-bend rigidity ratio to C/A = 1 and the stretching 
modulus to K/ksTa 2 = 10 3 , which gives negligible bond 
length fluctuations. The self-crossing of the chain is en- 
tirely prevented by setting tLj/ksT = 10. For sufficient 
numerical accuracy we choose a time step A = 0.0004. 
Output values are calculated every 10 3 -10 4 steps, total 
simulation times are 10 7 steps, the first 10 6 steps are not 
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included in the data analysis. The boundary condition 
at the forced end, d s r(0) = z, is realized by fixing the 
first two monomers in space by applying virtual forces, 
which also act (via the mobility tensor) on the rest of 
the filament. The rotational driving at the base imposes 
dt4>i = while force- and torque- free boundary con- 
ditions are adopted for the other end. The number of 
beads is in the range L/a = N = 30 — 100. Throughout 
this study, the persistence length is set to L p = 10a, thus 
L/Lp = 3-10. 

Buckling and plectonemes. To set the stage, 
we first repeat the zero-temperature scaling argument for 
the critical frequency lu c at which an axially rotated rod 
exhibits the buckling instability. At low rotational fre- 
quency, ujq < uj c , the rod is twisted but remains straight 
and the torque at the base balances the total rotational 
drag, ir-qa 2 uioL ~ Cf2(0). On the scaling level, the rod 
buckles when the twisting torque, Cf2(0), becomes com- 
parable to the bending torque, A/L, giving the critical fre- 
quency uj c ~ A I {-KT]a 2 L 2 ) independent of the twist rigidity 
C [15]. The linear stability analysis predicts an instabil- 
ity at tu c = 8.9A/(irr]a 2 L 2 ) [15], the numerical prefac- 
tor was confirmed by simulations in the zero-temperature- 
limit, i.e. for a very stiff polymer with persistence length 
L p = A/kgT L [16]. In Fig. 2 we show typical chain 
snapshots obtained in our dynamic simulations for driv- 
ing frequencies ujq/lu c — 0.3, 15 and 40 for a rather flexi- 
ble chain with L/L p = 10. For luq < lu c , the chain flexes 
randomly due to thermal motions and spins about its lo- 
cal axis at frequency ~ u>q. For luq ^> w c , in contrast, 
the polymer exhibits continuous generation of plectoneme- 
like structures and their diffusive transport from the driv- 
ing end to the free end. In the following we present a 
much-simplified scaling theory, valid for long elastic rods 
L ^> Lp, that establishes a minimal framework for treat- 
ing the competitive twist transport due to axial spinning 
and plectoneme creation. 

Scaling. — There are three different ways for the poly- 
mer to transport the injected twist to its free end, see 
fig. 1. The first one is the (A) axial spinning mode, where 
the polymer rotates around its contour like a speedometer 
cable, which we show to be the dominant dissipation mode 
for low enough driving frequencies. The second one is the 
(B) solid-body rotation mode, where the whole polymer 
coil whirls around the rotational axis at some frequency 
X- The third one is the (C) plectoneme creation/diffusion 
mode, in which plectoneme-like structures are continu- 
ously generated at the rotated end and are transported 
diffusively towards the free end. As we show in this paper, 
this dynamic twist-writhe conversion is a highly nonlinear 
mechanism that provides a very efficient means of reliev- 
ing torsional stress at elevated driving frequencies. Noting 
that in a stationary state, the twist that is injected into 
the polymer, ojq, has to exit the chain at the free end ei- 
ther in the form of axial spinning, solid-body rotation of 




Fig. 2: Typical chain conformations from our hydrodynamic 
simulations for chain length L — 100a and persistence length 
L p — 10a and rotational frequency ujq/uj c = 0.3, 15 and 40. 

frequency x, or writhe, we write 

lu = lu l + x + Aw, (2) 

where ujl is the axial spinning frequency of the free end 
and Auj denotes the fraction that is converted into writhe. 
In order to decide which of these three channels is in fact 
realized, we use the concept of minimum entropy produc- 
tion for non-equilibrium stationary states [20], according 
to which the state of least dissipation is stable. We will 
later confirm each of our scaling results by our simulations, 
which gives further credibility to our scaling approach. We 
thus have to estimate the power dissipation in each of the 
modes depicted in fig. 1. The power dissipation due to 
axial spinning is P as ~ r)La co^. Here we assume that av- 
erage axial spinning frequency is oil, meaning that in the 
plectoneme regime the rotational profile u>(s) decays very 
quickly (in fact exponentially) along the chain contour to 
the value ~ lu^. This is confirmed by our simulations, see 
fig. 4 (b). Likewise, the power dissipation due to solid- 
body rotation is given as P s b ~ r]R 3 x 2 - Except for a com- 
pact globule with R 3 ~ a 2 L we see that axial spinning is 
a less costly channel for twist transport than solid-body 
rotation, i.e. P as <C P s b- We therefore neglect solid-body 
rotation in what follows. The plectoneme-creation chan- 
nel is more complicated. We consider the creation of a 
single loop of radius R nearby the driven end, see fig. 1 
(C). Extension to a more complex plectoneme structure 
involving multi-loops (as actually observed in the simula- 
tions, see fig. 2) is straightforward, but does not change 
the conclusions on the scaling level. The bending energy of 
one loop is ksTLp/ R, the loop creation frequency is Au>, 
thus the power consumption is Pi oop ~ ksT 'Au)L p / R. In 
order to generate a loop of radius R nearby the driven 
end, the whole chain has to slide by the excess length 
2ttR during the time scale of (Auj)^ 1 . The frictional 
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Fig. 3: (a) Relative solid-body rotation frequency, x/luq, as a 
function of u>o/uj c . (b) Rescaled components of the vector Xi 
Xx = X± a ' 2 /(t J 'okBT), and Xz = X^o 2 l{lMsk B T), as a function 
of ujo/oJ a . Note that x = (xl + Xz) • The data are obtained 
for L = 50a and L„ = 10a. 



force associated with this sliding with velocity ~ RAuj 
is F s iide ~ r/LRAu;/ \n(L/a). Neglecting the logarithmic 
hydrodynamic correction, the power dissipation of sliding 
is Psiide ~ FshdeRAui ~ ijL(RAui) 2 . Minimizing the to- 
tal plectoneme creation dissipation, P p / ec = P( 00 p + Pslidei 
with respect to the loop radius i?, we obtain 



Pn, 



(k B TL p Aoj) A/3 (riL) 



a/3 



(3) 



with the loop size given by i? 3 ~ ksTLp/ (rjLAui). Due to 
the fractional power law P p i ec ~ (Aw) 4 / 3 , it is easy to see 
that plectoneme creation is unfavorable at low frequencies 
but will win over axial spinning at high frequencies. We 
now minimize the total dissipation P ~ P as + P p i ec with 
respect to the unknown frequency at the free end, w^, and 
use eq. (2) and X = 0. For low frequency we obtain pure 
axial spinning, lul ~ uq, i-e. the twist that is injected 
at one end comes out as axial spinning at the other end. 
For high frequencies coo > with a crossover frequency 
defined as 



— k B TL p /(Trr]a 3 L) ~ uj c (L/a), 



on the other hand, we obtain 



P 



2/3 1/3 
10*' UJ ' , 



rjLa u> 



2, 2/3, 4/3 



(4) 

(5) 
(6) 



Interestingly, the cross-over frequency u* is much larger 
than the zero-temperature twirling-whirling critical fre- 
quency uj c by a factor proportional to the chain contour 
length, L/a, and the dissipation P is reduced compared 
to the axial spinning scenario P as - 

To confirm these predictions, we extract the solid-body 
rotation rate vector, \, from our hydrodynamic simula- 
tions. For rigid-body rotation, i.e, dtY ~ \ x r, we ob- 
tain the vector x = (xx,Xz) via x = I 1 ' L, where 
I = J2jLii r jl ~ r j r j] i s the moment of inertia tensor and 
L = r x dr/dt is the polymer angular momentum 

(the mass of the polymer beads is set to unity) . Fig. 3b) 
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Fig. 4: (a) Steady-state profile of the rescaled twist density, 
fi(s) = fi(s)a, for two different driving frequencies, uio/u, — 
0.089 (squares) and 7.12 (triangles). The broken lines are the 
prediction from linearized theory, see text, (b) Steady-state 
profile of the rescaled rotational velocity, lo(s)/ujq together with 
an exponential fit (broken line), for a>o/w* = 3.56. All data 



were obtained for L = 50a and L v 



10a. 



shows x = (X-L)Xz) as a function of the rescaled driving 
frequency, loq/lu c ; solid-body rotation is small for cuq < lu c , 
while Xz grows significantly beyond u> c . The hydrody- 
namic shear due to this rotation deforms the plectonemes 
at high values of oj , which is a secondary effect that we 
neglect in the current version of our scaling arguments. 
In fig. 3 (a), we show that x/^o = |x|/ w o is much smaller 
than unity, which confirms that solid-body rotation is neg- 
ligible compared to other dissipation channels. 

The total twist Tw is defined as the integrated 
twist density along the chain arclength, Tw = 
(27r) _1 dsQ,(s). In the axial-spinning regime, i.e., coq < 
w*, chain shape fluctuations are decoupled from twisting 
motion, and the injected twist propagates diffusively along 
the chain [9]. The twist f2 thus obeys the linear diffu- 
sion equation, n r qa 2 dtVt = Cd 2 VL. In steady state under 
the boundary condition, CQ(0) = nr]a 2 u)oL, we thus ob- 
tain fl (s) = (irr]a 2 u>o/C)(s — L). This is in quantitative 
agreement with simulation for u>o/u>* = 0.089, as seen in 
fig. 4 (a) (open squares and broken line). The total twist 
in this regime is given by Tw ~ — irr/a 2 u) L 2 / (4irC) = 
— (A/ Ai:C)(L / a)(ujQ I uj*) and thus is linear in the driving 
frequency luq. In the plectoneme regime, on the other 
hand, the twist is much reduced compared to the linear 
law (lower broken line in fig. 4 (a)) and shows a nonlin- 
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ear spatial profile (open triangles in fig. 4 (a)). In this 
regime, the twist density O receives an additional con- 
tribution from the geometric torsion due to plectoneme 
formation, which however we neglect. Assuming the twist 
to be proportional to the average spinning frequency ojl 
in the region where plcctoncmes have already formed, i.e. 
Tw ~ (AL/aC)(ujL/^*) ~ (wq/w*) 1 / 3 , we thus predict 



' AL lu 
AiraC W* 



Tw 



(u>o < w») 



1/3 



(7) 



aC \ w, 



Finally, the effective rotational friction coefficient T r is 
defined via P ~ r r WQ, where P is the total dissipation, 
given by eq. (6) in the plectoneme regime. We obtain 

1 (w < W*) 

-2/3 (8) 



irr/a L 



For loq < w* , the rotational friction corresponds to the ax- 
ial spinning scenario, while for loq > w* a drastic decrease 
is predicted. Note that for ujq > cu* the dependence on 
chain length is T r ~ w 2 ^ 3 i 1 / 3 and thus increases only 
sub-linearly with polymer length. 

In Fig. 5, lul, Tw and T r from simulations are plotted 
as a function of ujq / and confirm the analytical predic- 
tions for the non-linear plectoneme regime , eqs. (5), (7) 
and (8). In the simulations, the torque Mq at a given 
driving frequency ujq is measured and averaged to obtain 
r r = (Mq)/u>q. In particular, the crossover frequency 
between the axial-spinning and the plectoneme regime is 
quite consistently found to occur at luq = c*w„ ~ (L/a)u> c 
with c* = 0.4 ± 0.1 and w* defined in eq. (4). Note that 
the whirling instability, at luq ~ lo c and realized only for 
L/Lp < 1, is conceptually distinct from the plectoneme 
transition at ujq ~ co*, which is only observable for long 
chains L/L p 3> 1. In fact, the two transitions do not 
merge or interconnect at intermediate values of L/L p : for 
L/L p ~ 1 a scmiflcxible chain rather shows a continuous 
shape and rotational mode evolution with increasing ujq 
without a sharply defined transition. 



Discussion. 



For a hydrodynamic diameter of ds- 



DNA of a ~ 2 nm, a bend persistence length L r 



30 



nm and length L ~ 12/im, we obtain according to eq. (4) 
a cross over frequency c^w* ~ 2 x 10 5 rad/s. The rota- 
tional friction of DNA molecules of length L ~ 12 /jm has 
been measured for rotational frequencies up to ujq ~ 12000 
rad/s in DNA unzipping experiments [10] and showed no 
detectable non-linear frequency dependence, in agreement 
with our estimate for the threshold On the other hand, 
for longer chains or in a crowded cellular environment with 
a much elevated viscosity, the crossover frequency a>* can 
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Fig. 5: Scaling plots of hydrodynamic simulation results, (a) 
Rescaled axial spinning frequency at free end, u>l/(*>*, (b) total 
twist divided by the chain length, \Tw\(a/L)), and (c) rescaled 
rotational friction constant, Y r / (ivqa 2 L) , plotted as a function 
of wo /w» . Soild lines denote the predictions from linear theory 
in the axial spinning regime valid for wo < w*, broken lines 
denote the non-linear results valid in the plectoneme regime 
as given in eqs. (5), (7) and (8). The crossover between the 
regimes occurs around a value of wo = c*w* with c* = 0.4±0.1. 



be lowered to an experimentally reachable value. The crit- 
ical torque M* = nr)a?Lc*u>* ~ c*7rfcsTL p /a is indepen- 
dent of both viscosity and chain length and for bare DNA 
is of the order of a few tens of ksT, too high for, e.g., single 
Escherichia coli RNA polymerase to achieve. But for chro- 
matin structures, the effective bending persistence length 
j^chrom nag ^ een ghown to stay rather constant while the 
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effective radius a chrom is increased substantially [6] , mean- 
ing that the critical torque M* om for a chromatin fiber 
might be of the order of ksT, quite in reach of typical 
torques generated by polymerase. 

In recent experiments, the elongation dynamics of su- 
pcrcoilcd DNA in response to sudden increase of tension 
was studied [14]. The dynamics showed almost no effect of 
the continuous removal of plcctoncmes, setting an upper 
bound on the hydrodynamic friction associated with DNA 
rotation, which turns out to be in agreement with the axial 
spinning scenario and thus is consistent with our results. 
Previous in vitro and in vivo studies on transcriptionally- 
drivcn DNA supercoiling have been interpreted as evi- 
dence of enhanced rotational friction when compared to 
the simple axial-spinning estimate [8]. Krebs and Dun- 
away showed in vivo that polymerase drives DNA super- 
coiling during the transcription process for linear DNA 
templates longer than 17 to 19 kbp [11]. One explanation 
involves static DNA bends or kinks stabilized by DNA- 
binding proteins, which are suggested to significantly slow 
down the torsional relaxation and lead to a considerably 
enhanced rotational drag torque [9, 12, 13]. It remains to 
study how static DNA bends interfere with the plectoneme 
creation/diffusion scenario developed in this paper. 

To summarize, using scaling arguments and hydro- 
dynamic simulations, we have studied the stationary 
noncquilibrium dynamics of torsionally driven scmiflcxiblc 
polymers. For chains much longer than their persistence 
length, L > L p , two dynamical regimes are distinguished: 
for small driving frequency, uj < , we find the standard 
axial-spinning regime where the chain flexes randomly and 
spins about its local axis. For large rotational frequency, 
ujq > j twist is locally converted into writhe close to the 
driven end and then diffuses out to the free end without 
much concerted solid-body rotation. In this plectoneme 
regime the filament exhibits only minimal axial spinning 
and a significant reduction of the rotational friction as 
compared to axial spinning is obtained. Two further con- 
clusions might be biologically relevant: The nature of the 
twist-writhe conversion process leads to a narrow spatial 
localization of the twist density close to the rotated part 
of the chain, which in-vivo might guide and concentrate 
the activity of twist-sensitive proteins to a region close to 
the polymerase complex. Finally, the positive supercoil- 
ing, created in front of the point of twist-injection in the 
form of plcctoncmes, would be a simple physical mecha- 
nism for loosening or even driving off histones from their 
nucleosomale core particles, the negative supercoiling be- 
hind the twist-injection point would form a template for 
strengthening or reforming the nuclcosomal core structure. 



REFERENCES 

[1] VOLOGODSKll A, Macromolecules, 27 (1994) 5623; 
Marko J. F., and SlGGlA E. D., Macromolecules, 28 

(1995) 8759; Hallatschek O., Frey E. and Kroy K., 
Phys. Rev. E, 75 (2007) 031905. 

[2] Ritort F, J. Phys.: Condens. Matter., 18 (2006) R531 

and references therein. 
[3] Marko J. F. and Siggia E. D., Phys. Rev. E, 52 (1995) 

2912. 

[4] Strick T. R., Allemand J.-F., Bensimon D., Bensi- 
mon A. and Croquette V., Science, 271 (1996) 1835. 

[5] Marko J. F., Phys. Rev. E, 55 (1997) 1758; Vologod- 
skii A. V. and Marko J.F., Biophys. J., 73 (1997) 123; 
Moroz J. D. and Nelson P., Proc. Nat. Acad. Sci. USA, 
94 (1997) 14418. 

[6] Lavelle C, Biochem. Cell. Biol., 87 (2009) 307. 

[7] Levinthal C. and Crane H., Proc. Natl. Acad. Sci. 
U.S.A., 42 (1956) 436. 

[8] Liu L. F. and Wang J. D., Proc. Natl. Acad. Sci. U.S.A., 
84 (1987) 7024; TSAO Y.-P., Wu H.-Y. and Liu L. F., 
Cell, 56 (1989) 111; Droge P. and Nordheim A., Nu- 
cleic Acids Research, 19 (1991) 2941. 

[9] Nelson P., Proc. Natl. Acad. Sci. U.S.A., 96 (1999) 
14342. 

[10] Thomen P., Bockelmann U. and Heslot F., Phys. 

Rev. Lett., 88 (2002) 248102; Nelson P., Phys. Rev. 

Lett, 92 (2004) 159801; Thomen P. and Heslot F., 

Phys. Rev. Lett, 92 (2004) 159802. 
[11] Krebs J. E. and Dunaway M., Mol. Cell. Biol., 16 

(1996) 5821. 

[12] Leng F. and McMacken R., Proc. Natl. Acad. Sci. 
U.S.A., 99 (2002) 9139; Leng F., Amado L. and Mc- 
Macken R., J. Biol. Chem., 279 (2004) 47564. 

[13] Stupina V. A. and Wang J. C, Proc. Natl. Acad. Sci. 
U.S. A, 101 (2004) 8608. 

[14] Crut A., Koster D. A., Seidel R., Wiggins C. H. 
and Dekker N. H., Proc. Natl. Acad. Sci. U.S.A., 104 
(2007) 11957. 

[15] Wolgemuth C. W., Powers T. R. and Goldstein R. 

E., Phys. Rev. Lett, 84 (2000) 1623. 
[16] Wada H. and Netz R. R., Europhys. Lett, 75 (2006) 

645. 

[17] Kamien R. D., Eur. Phys. J. B, 1 (1998) 1; Goldstein 
R. E., Powers T. R. and Wiggins C. H., Phys. Rev. 
Lett, 80 (1998) 5232. 

[18] Chirico G. and Langowski J., Biopolymers, 34 (1994) 
415. 

[19] Ermak D. L. and McCammon J. A., J. Chem. Phys., 
69 (1978) 1352. 

[20] DE Groot S. R. and Mazur P., N on- equilibrium Ther- 
modynamics (North-Holland Publishing Co., Amsterdam) 
1962, see Chap. V. 



* * * 



We thank C. Lavelle and P. Nelson for useful com- 
ments. Financial support from MEXT- Japan (Grant 
in Aid, No. 20740241) and the Excellence Cluster Nano- 
Initiative-Munich is acknowledged. 



